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O Abstract 

We present results for the universal anomalous dimension ~f U ni{j) of Wilson 
twist- 2 operators in the M = 4 Supersymmetric Yang-Mills theory in the first four 
orders of perturbation theory. 
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1 Introduction 



This paper deals with the study of the properties of the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [1] and Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [2] equations in the 
<-j \ JV = 4 Supersymmetric Yang-Mills (SYM) model [3]. 

Lev Nikolaevich Lipatov made a fundamental contributions in discovery of both the 
equations as in the framework of Quantum Chromodynamics (QCD) and later to its 
supersymmetric extensions. 
£S) | The BFKL and DGLAP equations resum, respectively, the most important contri- 

butions ~ a s ln(l/x) and ~ a s ln(Q 2 /A 2 ) in different kinematical regions of the Bjorken 
variable x and the "mass" Q 2 of the virtual photon in the deep inelastic scattering (DIS) 
and, thus, they are the cornerstone in analyses of the experimental data from lepton- 
nucleon and nucleon-nucleon scattering processes. 

In the supersymmetric cases the equations are simplified drastically and in the J\f = A 
SYM they become to be related each others for the nonphysical values of Mellin moments 
/\ • j as it has been proposed Lipatov in |4j. 

The purpose of this paper is to show the results for the anomalous dimension matrix 
of the twist- 2 Wilson operators, which have been obtained by author in collaboration with 
Lev Nikolaevich Lipatov during last 10 years. 

The anomalous dimensions govern the Bjorken scaling violation for parton distribu- 
tions in a framework of QCD. These quantities are given by the Mellin transformation 
(the symbol" is used for spin- dependent case and a s = a a /(47r)) 

la b {j) = C dx x^W^ a (x) = 7 W( J )a s + 7 «( J >2 + 7 ( 2 fe )( J ) a 3 + 0(a 4 )) 
Jo 

lab{j) = f'dxx^W^aix) = ^(j)a s + ^(j)a 2 s +^(j)a 3 s + 0(ai) (1) 

Jo 

of the splitting kernels W / 6_ >a (x) and Wb^ a (x) for the DGLAP equation [2j which evolves 
the parton densities f a (x,Q 2 ) and f a (x,Q 2 ) (hereafter a = A, g, (j) for the spinor, vector 



and scalar particles, respectively u) as follows 

alnLJ Jx y b 

7^fa(x,Q 2 ) = j 1 -Y,W^{x/y)Uy,Q 2 )- (2) 

dinQ 2 Jx y ^ 

The anomalous dimensions and splitting kernels in QCD are known up to the next-to- 
next-to-leading order (NNLO) of the perturbation theory (see [5] and references therein). 

The QCD expressions for anomalous dimensions can be transformed to the case of the 
TV-extended Supersymmetric Yang-Mills theories (SYM) if one will use for the Casimir 
operators CA,Cp,Tf the following values Ca = Cp = N c , TfUj = J\fN c /2. For J\f—2 
and M = 4-extended SYM the anomalous dimensions of the Wilson operators get also 
additional contributions coming from scalar particles [4]. These anomalous dimensions 
were calculated in the next-to-leading order (NLO) [U |6] for the M = 4 SYM. 

However, it turns out, that the expressions for eigenvalues of the anomalous dimension 
matrix in the M = 4 SYM [3] can be derived directly from the QCD anomalous dimensions 
without tedious calculations by using a number of plausible arguments. The method 
elaborated in Ref. [4] for this purpose is based on special properties of the integral kernel 
for the BFKL equation [H [TJ, [8] in this model and a new relation between the BFKL and 
DGLAP equations (see [4]). In the NLO approximation this method gives the correct 
results for anomalous dimensions eigenvalues, which was checked by direct calculations in 
Ref. [B|. Using the results for the NNLO corrections to anomalous dimensions in QCD [5] 
and the method of Ref. [4] we derive the eigenvalues of the anomalous dimension matrix 
for the J\f = 4 SYM in the NNLO approximation [9]. 

Starting from four loops, i.e. above existing QCD calculations, the corresponding 
results for the anomalous dimensions can be obtained (see [TQJ, [TTj, [T2]) from the long- 
range asymptotic Bethe equations together with some additional terms, so-called wrapping 
corrections, coming in agreement with Luscher approach. □ 

The obtained result is very important for the verification of the various assumptions 
(see recent reviews [H] p3] and references therein) coming from the investigations of the 
properties of a conformal operators in the context of AdS/CFT correspondence [T7] . 

The paper is organized as follows. In Section 2 we discuss the BFKL equation, the 
leading order anomalous dimensions of Wilson operators and propose the method of ob- 
taining the eigenvalues of the anomalous dimension matrix above the leading order. Sec- 
tion 3 contains the calculations of some Feynman diagrams by similar method. In Section 
4 we consider three-loop results for the universal anomalous dimension taking from the 
corresponding calculations in QCD. Four-loop corrections to the universal anomalous di- 
mension are considered in Section 5. 



■'■In the spin-dependent case a — A, g. 

2 The three- and four- loop results for the universal anomalous dimension have been reproduced (see 
[13] ) also by solution of so-called Baxter equation, which can be obtained from the long-range asymptotic 
Bethe equations. 



2 Evolution equation in J\f = 4 SYM 

The reason to investigate the BFKL and DGLAP equations in the case of supersymmetric 
theories is related to a common belief, that the high symmetry may significantly simplify 
their structure. Indeed, it was found in the leading order (LO) approximation [18], that 
the so-called quasi-partonic operators in M = 1 SYM are unified in supermultiplets with 
anomalous dimensions obtained from some universal anomalous dimension by shifting its 
argument by an integer number. Further, the anomalous dimension matrices for twist-2 
operators are fixed by the superconformal invariance [15] . Calculations in the maximally 
extended TV = 4 SYM, where the coupling constant is not renormalized, give even more 
remarkable results. Namely, it turns out, that here all twist-2 operators enter in the same 
multiplet, their anomalous dimension matrix is fixed completely by the super-conformal 
invariance and its universal anomalous dimension in LO is proportional to \l/(j — 1) — 
\I/(1) (see the subsection 2.2), which means, that the evolution equations for the matrix 
elements of quasi-partonic operators in the multicolour limit N c — y oo are equivalent to 
the Schrodinger equation for an integrable Heisenberg spin model [191 120] • In QCD the 
integrability remains only in a small sector of these operators [21] (see also [22J). In the 
case of M = 4 SYM the equations for other sets of operators are also integrable [231 121] • 
Similar results related to the integrability of the multi-colour QCD were obtained 
earlier in the Regge limit [25J. Moreover, it was shown [4], that in the M = 4 SYM 
there is a deep relation between the BFKL and DGLAP evolution equations. Namely, the 
j-plane singularities of anomalous dimensions of the Wilson twist-2 operators in this case 
can be obtained from the eigenvalues of the BFKL kernel by their analytic continuation. 
The NLO calculations in M = 4 SYM demonstrated [4], that some of these relations are 
valid also in higher orders of perturbation theory. In particular, the BFKL equation has 
the property of the hermitian separability, the linear combinations of the multiplicatively 
renormalized operators do not depend on the coupling constant, the eigenvalues of the 
anomalous dimension matrix are expressed in terms of the universal function j un i(J) which 
can be obtained also from the BFKL equation [4]. 

2.1 BFKL 

To begin with, we review shortly the results of Refs. [S], where the QCD radiative 
corrections to the BFKL integral kernel at t = were calculated, g We discuss only the 
formulae important for our analysis. 

The total cross-section a(s) for the high energy scattering of colourless particles A, B 
written in terms of their impact factors $j((&) and the t-channel partial wave G u (q, q') for 
the gluon-gluon scattering is 

. s f d 2 qd 2 q' , f a + l °° du /s\ w _ . ,. . .. ,. fo . 

"M- /(SOW * 4q) • a( « ) jL,SriW Gu{q ' q)t So = Mk| - (3) 

Here q and q' are transverse momentao of virtual gluons and s = 2paPb is the squared 
invariant mass for the colliding particle momenta pa and Pb- 



3 The t t^ case can be found in the recent papers [26] . 

4 To simplify equations hereafter we omit arrows in the notation of transverse momenta 



Using the dimensional regularization in the MS'-scheme to remove ultraviolet and 
infrared divergences in intermediate expressions, the BFKL equation for Gu{q,q*) can be 
written in the following form 



U)G u (q, qi ) = 6 D - 2 (q- qi ) + J d D - 2 q 2 K(q,q 2 )G ul (q 2} q 1 ] 



(4) 



where 

K(q u q 2 ) = 2u(q 1 )5 D - 2 (q 1 -q 2 ) + K r (q u q 2 ) (5) 

and the space-time dimension D = 4 — 2e for e — > 0. The gluon Regge trajectory oj(q) and 
the integral kernel K r (qi,q 2 ) related to the real particle production have been calculated 
in [27]-[29]. 

As it was shown in [7J [8], a complete and orthogonal set of eigenfunctions of the 
homogeneous BFKL equation in LO is 



2 \ 7-i 

y \ And 



Gn,j(q/q,0) = -« e 



<1 



(6) 



The BFKL kernel in this representation is diagonalized up to the effects related with 
the running coupling constant a s (q 2 ): 



«£" = 4a,tf 



AIS 



X(n, 7) + 5^L D (n, i)a s (q 2 ) 
Applying formulae of [8], we obtain the following results for eigenvalues (23): 
X (n,j) = 2*(l)-*( 7 + -)-*(l- 7 +-) (I 

^,7) = (y-2C(2)-^|)x(n ) 7) + 6C(3)H-^( 7 + f)+^(l-7+|) 

/ll 2 nA 1 9/ 
- 2$(n, 7 ) - 2<&(n, 1 - 7 ) - (^- - - -^J - X 2 (n, 7 ) 



(7) 



+ 



7T 2 COs(tT 7 ) 

sin 2 (7T 7 )(l — 2 7 ) 



(l + ^L) 7(1 ~ 7) .S 2 



, n,N 2 + 37(1-7) \ 5 o' 

" !L + iV3j(3-27)(l + 27)J \ 



(9) 



where <5™ is the Kroneker symbol, and ^(z), ty'(z) and ^"(z) are the Euler \1/ -function 
and its derivatives. The function $(n, 7) is given below 



\fc+i 



$(n ;7 ) = V^-J 

+ {-l) k ((3'{k + n + l) + (3'{k + l))^ 



q?'(k + n + l) -*'(*; + 1 

1 



fc + 7 + n/2 



*(A; + n + l) -*(jfe + i; 



(10) 



if, q , qt, q$, ... , i.e. in our formulae the momenta if, q , qt, q$, ... will be represented as 
q, q' , gi, <72, •■• , respectively. Note, however, that the momenta pa and ps are -D-space momenta. 



and 



?(*) = \ 



*>( Z -±±)-*>( Z ~ 



2 i ^2- 

Adding contributions of scalars and transforming fermions from fundamental to adjoint 
representation, we can obtain the BFKL form (JTj) in M = 4 SYM in DR scheme [30] 

<SV,7) = 6C(3) + *"( 7 + |)+*"(l- 7 + ^) 

2$(n,7)-2$(n,l- 7 )-2C(2)x(n,7), (11) 



where the DR coupling constant a s is related [3T] with the MS one a s as 

a s = a s + -a 2 s . (12) 

Note that the sum $(n, 7) + $(n, 1 — 7) can be rewritten (see [1]) as a combination 
of functions with argument dependent on 7 + n/2 = M and 1 — 7 + n/2 = M. Indeed 

$(n, 7) + $(n, 1 - 7) = X(n, 7) (£'(M ) + /?'(! - M)) 

+$ 2 (M) - /3'(M) [tf (1) - *(M)] + $ 2 (1 - M) - /3'(1 - M) [tf (1) - tf(l - M) 

where x( n >7) is given by Eq.© and 



* 2(M) = £ ftm £ FW ' (13) 

So, this transformation leads to the hermitian separability of BFKL equation in M = 4 
SYM (see Ref. [I] and discussions therein). 

2.2 Leading order anomalous dimension matrix in J\f = 4 SYM 

In the Af = 4 SYM theory [3] one can introduce the following colour and SU(4) singlet 
local Wilson twist-2 operators [U |6]: 

z SG^V^V^.-.V^^G^. , (14) 

z SG^V^V^.-.V^^G^. , (15) 

= S\? llMl V^...V H \ ai , (16) 

= 5^76 7w I>»...2? w A B< ) (17) 

= SftV^V^.V^r, (18) 

where T> p are covariant derivatives. The spinors \ and field tensor G PI1 describe gluinos 
and gluons, respectively, and r are the complex scalar fields. For all operators in 
Eqs. (I14p - (ll8p the symmetrization of the tensors in the Lorentz indices /ii,...,/ij and 
a subtraction of their traces is assumed. 



Q9 


.,flj 


Q9 


.,Hj 




.,Hj 




.,Hj 


Q<t> 


.,Hj 



The elements of the LO anomalous dimension matrix in the M = 4 SYM have the 
following form (see [20J): 

for tensor twist-2 operators 



iSu) 



4 [ *(1) - *(j - 1) - - + 



ifU) 



1 



2 2 

+ 






J 3 + 1 J + 2 / 

2 2 1 

- + - 



J 3 + 1 3 + 1, 
lf 9 {3) - 



12 



J + l j + 2. 



\3-l 3 3 + 1 J 



3 3 + 1, 



7S0') = 4(M/(1) -M>( J + l)) 






7$(j) 



J + l' 



J ~ 1 3. 



(19) 



for the pseudo-tensor operators: 



«(i) 



4 (*(1) - *(7 + 1) " J^ + J 



tSj) 



1 2 

+ 



V 3 3+1 J 
4 U(l) - *(j + 1) + 



7?(i) 



J' J + 1 . 



J + l J 



(20) 



The matrices, based on the anomalous dimensions ( fl9l) and ( !20l) . can be diagonalized 
20, 4j. They have the following remarkable form 



LTZT 1 



Af=4 



unpol 



-45 x 0'-2) 

-4-5x0-) 

-ASxU + 2) 



DTD 



-i 



iV=4 



- pol 



o -4^0 + 1; 



where Si(J) is defined below in (|25|) . 

Thus, the LO anomalous dimensions of all multiplicatively renormalized operators can 
be extracted through one universal function 



7 (0) (?) 



AS(j - 2) = -4$(j - 1) - *(1)) = -4 £ 



i-2 



r=l 



2.3 Method of obtaining the eigenvalues of the anomalous di- 
mension matrix in M = 4 SYM 

As it was already pointed out in the Introduction, the universal anomalous dimension 
can be extracted directly from the QCD results without finding the scalar particle con- 
tribution. This possibility is based on the deep relation between the DGLAP and BFKL 
dynamics in the J\f = 4 SYM 0]. 

To begin with, the eigenvalues of the BFKL kernel are the analytic functions of the 
conformal spin \n\ at least in two first orders of perturbation theory (see Eqs. (J7J), (|SJ) 
and (ITT]) ). Further, in the framework of the DR-scheme [30] one can obtain from (JSJ) and 
(jSD, that there is no mixing among the special functions of different transcendentality 
levels i lj, i.e. all special functions at the NLO correction contain only sums of the terms 
~ 1/Y {i = 3). More precisely, if we introduce the transcendentality level % for the 
eigenvalues a; (7) of integral kernels of the BFKL equations in an accordance with the 
complexity of the terms in the corresponding sums 

tf ~ I/7, *' ~ 0' ~ C(2) ~ 1/7 2 , *" ~ 0" ~ $ ~ C(3) ~ 1/7 3 , 

then for the BFKL kernel in LO and in NLO the corresponding levels are i — 1 and i — 3, 
respectively. 

Because in M = 4 SYM there is a relation between the BFKL and DGLAP equations 
(see [HI H]), the similar properties should be valid for the anomalous dimensions them- 
selves, i.e. the basic functions 7„ ni (j), luniii) an d llniU) are assumed to be of the types 
~ \j f with the levels i — 1, i — 3 and i — 5, respectively. An exception could be for the 
terms appearing at a given order from previous orders of the perturbation theory. Such 
contributions could be generated and/or removed by an approximate finite renormaliza- 
tion of the coupling constant. But these terms do not appear in the DR-scheme. 

It is known, that at the LO and NLO approximations (with the SUSY relation for 
the QCD color factors Cp = Ca = N c ) the most complicated contributions (with % — 1 
and 2 = 3, respectively) are the same for all LO and NLO anomalous dimensions in 
QCD |5J and for the LO and NLO scalar-scalar anomalous dimensions [6]. This property 
allows one to find the universal anomalous dimensions 7„„j(j) and 7„ n j(j) without knowing 
all elements of the anomalous dimensions matrix [4], which was verified by the exact 
calculations in [B]. 

Using above arguments, we conclude, that at the NNLO level there is only one possible 
candidate for 7„ n j(j). Namely, it is the most complicated part of the QCD anomalous 
dimensions matrix (with the SUSY relation for the QCD color factors Cp = Ca = N c ). In- 
deed, after the diagonalization of the anomalous dimensions matrix its eigenvalues should 
have this most complicated part as a common contribution because they differ each from 
others only by a shift of the argument and their differences are constructed from less com- 
plicated terms. The non-diagonal matrix elements of the anomalous dimensions matrix 
contain also only less complicated terms (see, for example, anomalous dimensions exact 
expressions at LO and NLO approximations in Refs. [5] for QCD and [6] for Af = 4 SYM) 



5 Similar arguments were used also in 32 to obtain analytic results for contributions of some compli- 
cated massive Feynman diagrams without direct calculations (see also the section 3). 



and therefore they cannot generate the most complicated contributions to the eigenvalues 
of anomalous dimensions matrix. 

Thus, the most complicated part of the NNLO QCD anomalous dimensions should 
coincide (up to color factors) with the universal anomalous dimension 7„ ni (j). 



3 Calculation of Feynman integrals 

Similar arguments give a possibility to calculate a large class of Feynman diagrams, so- 
called master- integrals [33]. Let us consider it in some details. 

Application of the integration-by-part (IBP) procedure [31] to loop internal momenta 
leads to relations between different Feynman integrals (FI) and, thus, to necessity to cal- 
culate only some of them, which in a sense, are independent (see [35] ). These independent 
diagrams (which were chosen quite arbitrary, of course) are called the master-integrals 
[33]. 

The application of the IBP procedure [31] to the master-integrals themselves leads to 
the differential equations [36 J37] for them with the inhomogeneous terms (ITs) containing 
less complicated diagrams. U The application of the IBP procedure to these diagrams 
leads to the new differential equations for them with the new ITs containing even farther 
less complicated diagrams. Repeating the procedure several times, at a last step one can 
obtain the ITs containing only tadpoles which can be calculated in-turn very easily. 

Solving the differential equations at this last step, one can reproduce the diagrams for 
ITs of the differential equations at the previous step. Repeating the procedure several 
times one can obtain the results for the initial Feynman diagram. 

This scheme has been used successfully for calculation of two-loop two-point [351 EH] 
and three-point diagrams [3_Hl [321 HO] with one nonzero mass. This procedure is very 
powerful but quite complicated. There are, however, some simplifications, which are 
based on the series representations of Feynman integrals. 

Indeed, the inverse-mass expansion of two-loop two-point and three-point diagrams [j 
with one nonzero mass, can be considered as 



FI 



Q 



k^i»n 



n=l 



rr 



ln(-x)F lil (n) + -F li2 (n) 



+ 
+ 



1 



\n\-x) F 2il (n) + - ln(-a;) F 2>2 (n) + - F 2>3 (n) + C(2) F 2A (n) 



\n-\-x) F 3jl (n) + - \n 2 (-x) F 3 , 2 (n) + - ln(-x) F 3 , 3 (n) + - F 3A (n) 



+C(2)ln(-x)F 3 , 5 H + C(3)F 3 , 6 (n) 



+ 



(21) 



where x = q 2 /m 2 , 77 = 1 or — 1, c = 0, 1 and 2, and a = 1 and 2 for two-point and 
three-point cases, respectively. 



6 The "less complicated diagrams" contain usually less number of propagators and sometimes they can 
be represented as diagrams with less number of loops and with some "effective masses" (see, for example, 
[38] and references therein). 

7 We consider only three-point diagrams with independent momenta q\ and (72, which obey the condi- 
tions q\ = q 2 



Oand ( qi +q 2 ) 2 =q 2 ^ 0. 



Here the normalization N = (~p 2 /m 2 ) £ , where /Z = 47re~ TB / u is in the standard MS- 
scheme and 7^ is the Euler constant. Moreover, the space-time dimension is D = 4 — 2e 
and 

C n = 1 (22) 

for diagrams with one-massive-particle-cuts (m-cuts) and 

C n = 1, and C n = ^- = C n (23) 

for diagrams with two- massive-particle-cuts (2m-cuts). 

For m-cut case, the coefficients F^^{n) should have the form 

F N , k (n) ~ %^. (24) 

In this section S± a = S± a (j - 1), S± a> ± b = S± a> ± b (j ~ X )' s ±a,±b,±c = S± a>±bt±c (j - 1) are 
harmonic sums 

3 (±i)m 3 (±i)m 

s± a (j) = E^> s ±at±b , ±c ,...(j) = E^-fVc-H, (25) 

m=l " L m=l " b 

For 2m-cut case, the coefficients .F/v,fc(n) should have the formal 

t? 1 \ ^±0,... Ki,... Vv a . / oc x 

i*V,fc(n) — , — r , — — 26 

n" n 6 n 6 



where 



W) = E% ^ c ,.(j) = E §^ C ,.M, (27) 



m=l m m=l m 



W.tt) =tS WW,...(;) = t ^S btC ,..(m), (28) 



m=l m mTl m 



The terms ~ V a ,... and ~ W a ,... can come only together with the coefficients C n — 1 
and C n = C* n , respectively. The terms ~ S± aj _ can appear in combination with both 
C n values. The origin of the appearance of the terms ~ V a ,... and ~ W 0j „. in the 2m-cut 
case, is the product of series (I2T]) with the different values of the coefficients C n = 1 and 

r - r 

As examples, consider two- loop two-point diagrams Ii, I 5 and ii 2 , studied in [32] 

Ii = ^E-U^ 2 (-^)--ln(-^)+C(2) + 25 2 -2^ + A}, (29) 

q ~l n 12 n n n z ) 

h = ^E^ Z ^(- ln2 (-^) + -M-^)-2C(2)-4^ 2 -^-2^j,(30) 
q ~l n I n n z n J 

h2 = ^E^{ 1 + 7^(-21n(-x)-3iy 1 + -)}. (31) 

<T ~^ n 2 in {2n)\ \ nj ) 



8 Really, there are even more complicated terms as ones in Eqs. (58) and (59) of [32], which come from 
other 77 values in (|21j) . However, they are outside of our present consideration. 

9 



From (|29|) and (!30|) one can see that the corresponding functions Fj^^{n) have the 
form 



F N ,k( n ) 



1 



n 



2-iV' 



(iV > 2), 



(32) 



if we introduce the following complexity of the sums (X)£Li a « = a ) 



$»?a ~ ^Tjai^aa ~ § r, ai , n a 2 ,■ ■ ■ ,qa m ~ Ca ~ — , 



(33) 



where $ = (S,V,W). 
In Eq. ([S]), 



F N A n ) ~ 



1 



/? 



1-iV' 



(iV>l) 



(34) 



since now the factor 1/n 2 has been already extracted. 

So, Eqs. (I29]) - (I3~T]) show that the functions F^^{n) should have the following form 



1 r, / X 1 

F N h{n) ~ 



ir 



n 



3-7V ' 



(N>2) 



(35) 



and the number 3 — N defines the level of transcendentality (or complexity) of the co- 
efficients FN,k{n). The property reduces strongly the number of the possible elements 
in FN,k(n). The level of transcendentality decreases if we consider the singular parts of 
diagrams and/or coefficients in front of (-functions and of logarithm powers. 

Other /-type integrals in [32] have similar form. They have been calculated exactly 
by differential equation method [36j [37] . 

Now we consider two-loop three-point diagrams, Pi, P 5 , P 6 , P i3 and Pi 2 , considered 
in 



Pi 



N 



(Q 



2\2 



2_^ m 1 o^3 ^2 """ O,- 



n=l 



n 



2e 3 e 2 



2s 



55 2 -S\ + - 2 -- ln(-x) + \ ln 2 (-x) - C(2) 



-3C3 



Si + - C2 + -S3 + -S1S2 + -5? + 4^ + 24 + — 
71/ 3 2 n n z n A 



C2 _45 2 _ 2 ^-A)b,(-,r) 
n rr 



S '+2^ ln 



-x 



^ln 3 (-x) 



iV 



(9 



2^2 



E 

n=l 



77 



-6C3 + 2(5xC 2 + 65 3 - 2S 1 S 2 + 4^-^ + 2^ 

n n rr 



+ (-45 2 + 5 2 - 2=±) ln(-x) + S ± ln 2 (-x)| 






E 



-x 



?? 



ln(— x) 



n 



C2 - 35 2 - 45_ 2 - 3 



Si 



n rtr 



+ f35 1 + -Vn(-x)-|ln 2 (- i 
\ nj 2 

+105_ 3 - 125.3,! - 45 x 5_ 2 - l^ - ?S - 5^ - 1 + f ~5 2 - !U 2 



f2C 3 4 (7S'i + - 1(2-253-95x52 



2 71 2 n ~n 2 n 3 ' V2 ' 2" 



(36) 



(37) 



10 



^ + - 2 - 2Q 2 ) ln(-x) + \ (l S x + -) ln 2 (-x) + \ ^ 
n n z / 2\ nj o 



+5— + — - 2C 2 ln(-x) + - [7S 1 + - ln 2 (-x) + - ln 3 (-a;) , (3 



iV ^ „ f S 2 1 



Pl3 (? 2 ) 2 £ X I 2e 2 2e 



S 3 + 45! 2 - 4^ 



+ yC 2 



-Si 



,3 



35-3,1 + 351,1,2 + 35-1,2,! - Si + (jS 3 - 8Si, 2 )Si + ^S 2 S 2 } , (39) 

N r n (nl) 2 r 2 2 / 1 \ 

p » - 7^m^ + -A s >- 3w i + -n- h ^) +i2w >- i8w " 

-135 2 + S 2 - 65!Wi + 2^ + A _ 2^1 + -") ln(-x) + ln 2 (-x)l , (40) 

Now the coefficients FM,k{n) have the form 

^,*(n) ~ -1^, (iV>3), (41) 

The diagrams Pi, P 5 and P 6 (and also P 3 in [32]) have been calculated exactly by 
differential equation method [35J[37]. 

To find the results for P 13 and P12 (and also all others in [32J) we have used the 
knowledge of the several n terms in the inverse-mass expansion (J2"T1) (usually less than 
n = 100) and the following arguments (see [40] and discussions therein): 



The coefficients should have the structure 04 ip with the rule (133]) . The condition 04 ip 
reduces strongly the number of possible harmonic sums. It should are related with 
the specific form of the differential equations for the considered master integrals, 
like 

( ke + m 2 - — - FI = less complicated diagrams , 
\ dm z ) 

with some k values. We note that for many other master integrals (for example, 
for sunsets with two massive lines in [351 E]) the property fj4T|) is violated: the 
coefficients Fj^^in) contain sums with different levels of complexity, u 

If a two-loop two-point diagram with the "similar topology" (for example, Ii for 
Pi and P3, J5 for P 5 and Pq, I\i for P i2 an so on) has been already calculated, we 
should consider a similar set of basic elements for corresponding F^^(n) of two-loop 
three-point diagrams but with the higher level of complexity. 

Let the considered diagram contain singularities and/or powers of logarithms. 

Because in front of the leading singularity, or the largest power of logarithm, or 
the largest (^-function the coefficients are very simple, they can be often predicted 
directly from the first several terms of expansion. 



9 Really, Refs. [35j [41] contain the Nilson polylogarithms, whose sum of indices relates directly to the 
level of transcendentality (4 — N). The representation of the series (f29|) -([3T j) and (j3"6]) -([4"0 |) . containing 
S± a ,..., as polylogarithms can be found in [32 for m-cut case and in |42j for 2m-cut one, respectively. 
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Moreover, often we can calculate the singular part using another technique (see 
for extraction of ~ W\(n) part). Then we should expand the singular parts, find 
the basic elements and try to use them (with the corresponding increase of the level 
of complexity) to predict the regular part of the diagram. If we have to find the 
£-suppressed terms, we should increase the level of complexity for the corresponding 
basic elements. 

Later, using the ansatz for F^^{n) and several terms (usually, less than 100) in the 
above expression, which can be calculated exactly, we obtain the system of algebraic equa- 
tions for the parameters of the ansatz. Solving the system, we can obtain the analytical 
results for FI without exact calculations. To check the results, it is needed only to calcu- 
late a few more terms in the above inverse-mass expansion ( l2Tj) and compare them with 
the predictions of our anzatz with the above fixed coefficients. 

The arguments give a possibility to find the results for many complicated two-loop 
three-point diagrams without direct calculations. Some variations of the procedure have 
been successfully used for calculating the Feynman diagrams for many processes 

(see [SHI E21 EH1 S33)- 

4 Universal anomalous dimension for J\f = 4 SYM 

The final three-loop result lj for the universal anomalous dimension j U ni{j) for Sf = 4 
SYM is [9] 

l(j) = luniij) = a,l$!b(j) + a 2 luliti) + a 3 luniU) + -, a = -^- , (42) 

~720' + 2) = -S u (43) 

~7inl(i + 2) = (S 3 + S. 3 )-2S. 2 , 1 + 2S 1 (S 2 + S. 2 ) 7 (44) 

^7S(j + 2) = 2 S- 3 S 2 - S 5 - 2 S„ 2 S 3 - 3 3L 5 + 24 3_ 2)liM 



where 



+6 [ S'_4 ) i + S_3 :2 + S-2,3 j ~ 12 [ S_3,l,l + '5-2, 1,2 + "5-2,2,1 

5 2 + 2 sf\ (3 S_ 3 + S 3 - 2 SL 2il ) - S t (8 S_ 4 + S 2 _ 2 

+4 5 2 5_ 2 + 2 S 2 2 + 3S 4 - 12 5-3,1 - 10 5_ 2 , 2 + 16 5_ 2 ,i,i) (45) 

and S a = S a (J), S a)b = S ayb (j), S aAc = S a ^ c (j) are harmonic sums (see Eq. ([25]) and 

5- aAc ,...(j) = (-l) j S- a , b , C7 ...(j) + S- a7b , c ,..(oo)(l- (-l) j ). (46) 

The expression ( 1461) is defined for all integer values of arguments (see [HI HJ H5] ) but can 
be easily analytically continued to real and complex j by the method of Refs. [HJ I4"5|l4"6]. 



10 Note, that in an accordance with Ref. [7 our normalization of j(j) contains the extra factor —1/2 
in comparison with the standard normalization (see [4]) and differs by sign in comparison with one from 
Ref. 0. 
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4.1 The limit j -> 1 

The limit j — > 1 is important for the investigation of the small-x behavior of parton 
distributions (see review [47] and references therein). Especially it became popular re- 
cently because there are new experimental data at small x produced by the HI and ZEUS 
collaborations in HERA |48j . 

Using asymptotic expressions for harmonic sums at j — 1 + u — > 1 (see [H |9]) we 
obtain for the M = 4 universal anomalous dimension j un i(j) in Eq. (T42J) 

7 £i(l+w) = -32C8 + 0(w 1 ), (48) 



7^(1 +w) = 32C 3 ^-232C4--112OC 5 + 256C3C2 + 0(w 1 ) (49) 

in an agreement with the predictions for 7^(1 + 0;), 7^(1 +w) and also for the first term 
of 7„ n j(l + u) coming from an investigation of BFKL equation at NLO accuracy in [Sj. 

4.2 The limit j -► 4 

The investigation of the integrability in A/" = 4 SYM for a BMN-operators [49] gives a 
possibility to find the anomalous dimension of a Konishi operator [501124"] . which has the 
anomalous dimension coinciding with our expression (1421 for j = 4 



7«m(j)li=4 = -6a s + 24 a 2 s - 168 af (50) 

It is confirmed also by direct calculation in two [HU E] and three-loop [52] orders. The 
four and five loop corrections to the anomalous dimension of a Konishi operator have 
been also calculated recently in [53], [54] and [55], [56], respectively (see the recent review 
[T6] and references therein). 



4.3 The limit j — )• 00 

In the limit j — > 00 the results (14"3"I) - (14"5l) are simplified significantly. Note, that this limit 
is related to the study of the asymptotics of structure functions and cross-sections at 
x — > 1 corresponding to the quasi-elastic kinematics of the deep-inelastic ep scattering. 

We obtain the following asymptotics for the M = 4 universal anomalous dimension 
luniU) in Ec l- (H2J) with 

1 { 1U) = -4(lnj+7i.)+0(j- 1 ), (51) 

7&0') = 8C 2 (lnj+ 7s ) + 12C3 + 0(j- 1 ), (52) 

720') = -88C 4 (lnj+ 7 i?)-16C2C3-80C5 + O(j- 1 ), (53) 

where 7^ is Euler constant (see also the normalization in Eq. ( 12TT) ). 
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4.3.1 Resummation of ^ uni and the AdS/CFT correspondence 

Last several years there was a great progress in the investigation of the M = 4 SYM 
theory in a framework of the AdS/CFT correspondence [17] where the strong-coupling 
limit a s — > oo is described by a classical supergravity in the anti-de Sitter space AdS 5 x S 5 . 
In particular, a very interesting prediction [57] (see also [58]) was obtained for the large- j 
behavior of the anomalous dimension for twist-2 operators 

aN c 

jU) = a ( z ) ln J » z = — - = 4a « ( 54 ) 

TV 

in the strong coupling regime (see Ref. [59] for asymptotic corrections): 

lim a = -z 1 ' 2 + — + O (V 1/2 ) . (55) 

z->oo 37J- V / y ' 

On the other hand, the results for ■y un i(j) in Eqs. (142]) and (|5T]) - (!53|) allow one to find 
three first terms of the small- z expansion of the coefficient a(z) 

71-2 11 

lima = -zH z 2 - -— ttV + ... . (56) 

*->o 12 720 v ; 

For resummation of this series Lipatov suggested the following equation for the ap- 
proximation a [B] 

71-2 

z = —a H a 2 (57) 

12 v ; 

interpolating between its weak-coupling expansion up to NNLO 

a = -z + — z 2 - — vr 4 z 3 + 0(z A ) (58) 

12 72 v ; v ; 

and strong-coupling asymptotics 

~ a = _^! Z W + ± + iz- 1 ' 2 ) « -1.1026 ^ 1/2 + 0.6079 + O (V 1/2 ) . (59) 

It is remarkable, that the predictions for NNLO based on the above simple equation and 
obtained before the NNLO results (Hoi) and ( [ 53 ] )). are valid with the accuracy ~ 10%. It 
means, that this extrapolation seems to be good for all values of z. I 11 ! 

4.3.2 Beisert-Eden-Staudacher equation 

Recently the integral Beisert-Eden-Staudacher (BES) equation has been proposed in [60] 
for some function f(x), which is related with a(z) of (|54|) at x = 0, i.e. /(0) = a(z). 
At small coupling constant z, this equation gives a lot of coefficients c m of the expansion 

/(o) =Ec m ^ m . 

m=0 



11 Some improvement of (J57I) van be found in [61] . 
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These coefficients c m obey to the transcendentality principe, i.e. c m ~ C(2w) for m > (or 
products of (- function with the sum of indices equal to 2m). Moreover, up to 4-loop, the 
coefficients are in agreement numerically with ones, obtained directly from calculations 
of Feynman diagrams [2H E2] • 

The most important purpose, however, is to find the z — > oo limit from the BES 
equation, i.e. to try to reproduce the Polyakov et al. asymptotics ~ z 1 ^ 2 (see the r.h.s. 
of fl55l) ). The study was performed and the asymptotics were reproduced numerically [63] 
and analytically [M] . 

Recently authors of [65] found a method to evaluate the c m coefficients of the expansion 



/(o) = E 



i^rn. Z 



(l-m)/2 



m=0 



of the BES equation and calculated several of them. The first three coefficients are in 
agreement with the results of exact calculations performed in [57], [5H] and [66], respec- 
tively. Moreover, the results of [65J are in well agreement with transcendentality principe: 
C\ ~ In 2 and c m ~ C( m ) for m > 1 (or products of (^-function with the sum of indices 
equal to m). 

5 Bethe-ansatz and four-loop universal anomalous di- 
mension 

The long-range asymptotic Bethe equations for twist-two operators have the form 



x 



M - _ + /-i _ 2/_+„-\ M H 

■^k ■L'm \ L if l^k^m) _ ^ in int. . \\ TT * 



n-^fc -^m V- 1 - tl ^k^m) / ■ at \\ TT k i tan\ 

w - t + -t- (l-ayr-r+) ex P( 2 ^KO), H^ = 1 -( 60 ) 

\ k / rn=l,m^k k ^m \ L i) I Jj k Jj m) k=\ k 

These are M equations for k — 1, . . . , M Bethe roots «&, which need to be solved for the 
Bethe roots «&. The variables x^ are related to Uk through Zhukovsky map 




xt = x(ut), u ± = u±-, x(u) = - M+Wl-4^ , (61) 



The dressing phase 6 ~ ((3) is a rather intricate function conjectured in ^U\ and its exact 
form is not so important for the present consideration. 

Once the M Bethe roots are determined from above equations for the state of interest, 
its asymptotic all-loop anomalous dimension is given by 

7 ^(,) = 2, 2 f;(^-^y (62) 

fc=i \ x k x k ) 

The above equations can be solved recursively order by order in g at arbitrary values of 
M once the one-loop solution for a given state is known. 

This technical problem can nevertheless be surmounted. Assuming the maximum tran- 
scendentality principle [I] at four-loop order one may derive the corresponding expression 
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for the anomalous dimension by making an appropriate ansatz with unknown coefficients 
multiplying the nested harmonic sums, and subsequently fixing these constants. The lat- 
ter is done by fitting to the exact anomalous dimensions for a sufficiently large list of 
specific values of M as calculated from the Bethe ansatz. I 12 l 

Luckily, at one-loop the exact solution of the Baxter equation is known [57] and is 
given by a Hahn polynomial. Knowing the one-loop roots one can then expand equation 
( 160]) in the coupling constant g order by order in perturbation theory. The equations for 
the quantum corrections to the one-loop roots are of course linear, and thus numerically 
solvable with high precision. 

The result for the four-loop asymptotic dimension has the form [TO] (M = j + 2): 

^aba {j + 2) = 4 :S _ 7 + 6S 7 + ... + -C{3)S 1 {S 8 -S- a + 2S- V ), (63) 

where the symbol ... marks large set of the nested sums of degree seven. 

It is possible to analytically continue the expression in the r.h.s. of (I6"3"j) to the vicinity 
of the pomeron pole at M = — 1 + u. An explanation for how this is done may be found 
in [15]. 

Harmonic sums of degree seven may lead to poles no higher than seventh order in u. 
In fact, it is known that none of the sums in r.h.s. of (|63l) can produce such a high-order 
pole except for the two sums 5*7 and S--?. Their residues at 1/u 7 are of opposite sign. 
Thus, one immediately sees that the sum of the two residues does not cancel. 

However, from BFKL calculations [H H] , it is possible to conclude that at the vicinity 
of the pomeron pole at M = — 1 + u the four loop anomalous dimensions 

7«m(l+w) ~ 1/cA (64) 

It proves, that the above result is not full and there are so-called wrapping corrections. 
The contribution of the wrapping corrections has been added in [11] . So, the full result 
has the following form 



7«n*0' + 2) 7^ A (j + 2)+7™(j + 2), 

2 5*_5 + 2 S5 + 4 (S^i — 03 : _2 + 5'_2 I _3 — 2 D_2,_2,l; 



1 7^(j + 2) = I Si 



256 ,umw ' 2 

-4S_ 2 C(3)-5C(5) 

This result is in full agreement with BFKL predictions ([61 

We note that using similar technique and a property of reciprocity (see [68J and ref- 
erences therein), the five-loop corrections fo universal anomalous dimensions have been 
found in 



6 Conclusion 

In this review we presented the anomalous dimension 7 nTl i(j) for the M = 4 supersym- 
metric gauge theory up to the next-to-next-to-next-to-leading approximation. All the 



12 The study is similar to one considered in the Section 3 and used for calculations of Feynman integrals. 
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results have been obtained with using of the trans cendentality principle. At the first three 
orders, the universal anomalous dimension have been extracted from the corresponding 
QCD calculations. The results for four and five loops have been obtained from the long- 
range asymptotic Bethe equations together with some additional terms, so-called wrapping 
corrections, coming in agreement with Luscher approach. 

A.V.K. thanks to Binur Shaikhatdenov for careful reading of the paper. 
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